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ABSTRACT 

Hot accretion tori around a compact object are known to be susceptible to a global 
HH ■ hydrodynamical instability, the so-called Papaloizou-Pringle (PP) instability, arising 

| from the interaction of non-axisymmetric waves across the corotation radius, where the 

O |. wave pattern speed matches the fluid rotation rate. However, accretion tori produced 

in various astrophysical situations (e.g., collapsars and neutron star binary mergers) 
£^ ■ are likely to be highly magnetized. We study the effect of magnetic fields on the PP 

c/3 ' instability in incompressible tori with various magnetic strengths and structures. In 

& . general, toroidal magnetic fields have significant effects on the PP instability: For thin 

tori (with the fractional width relative to the outer torus radius much less than unity), 
(f) . the instability is suppressed at large field strengths with the corresponding toroidal 

J> ' Alfven speed va<$> il 0.2rfl (where is the flow rotation rate). For thicker tori (with 

the fractional width of order 0.4 or larger), which are hydrodynamically stable, the 
instability sets in for sufficiently strong magnetic fields (with va</> <; 0.2rQ). Our results 
. suggest that highly magnetized accretion tori may be subjected to global instability 

even when it is stable against the usual magneto-rotational instability. 
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1 INTRODUCTION 



. Differentially rotating fluid flows, ubiquitous in astrophysics, can exhibit rich dynamical behaviors. Papaloizou & Pringle 
(1984) discovered that accretion tori can be subjected to a global non-axisymmetric instability that grows on a dynamical 
time-scale. Accretion tori are bagel-shaped discs with high internal temperatures and well-defined boundaries. They may be 
representative of certain stages or regions of the inner accretion flows around black holes, such as those found in active galactic 
nuclei and quasars (e.g., Begelman, Blandford & Rees 1984). They may also form in the gravitational collapse of the rotating 
core of massive stars (e.g., Woosley 1993) and after the merger of compact neutron star and black hole binaries (e.g., Duez 
et al. 2009; Etienne et al. 2009; Rezzolla et al. 2010; Montero et al. 2010), and thus are thought to be the central engine of 
gamma-ray bursts (e.g., Meszaros 2006). The Papaloizou-Pringle (PP) instability arises from the interaction between non- 
axisymmetric waves across the corotation radius (r c ), where the wave pattern rotation frequency f2 p equals the background 
fluid rotation rate (e.g., Blaes & Glatzel 1986; Goldreich, Goodman & Narayan 1986; Glatzel 1987b). Waves outside the 
corotation radius (r > r c ) have Q, p larger than f2(r) and carry positive energy, while waves at r < r c have Q p < fi(r) anc l carry 
negative energy. Instability occurs when the negative-energy waves inside r c lose energy to the positive-energy waves outside 
r c , leading the amplification of the wave amplitudes. To maintain the instability, the waves must be efficiently reflected at 
the inner and outer boundaries so that they are trapped in the torus. The growth rate of the PP instability is maximal for a 
constant-angular momentum torus. For a very thin torus (with the inner and outer radii close to each other), the instability 
disappears when p — dlnfi/dmr > — v3; for a wider torus, the instability persists (with decreasing growth rate) as the 
Keplerian rotation profile (p = —3/2) is approached (e.g., Papaloizou & Pringle 1985, 1987; Goldreich et al. 1986; Zurek 
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& Benz 1986; Sekiya & Miyama 1988). Other properties of the PP instability, such as its connection with the instability of 
vortices (e.g., Glatzel 1987a), its non-linear evolution (e.g., Goodman, Narayan & Goldreich 1987; Hawley 1991) and the effect 
of accretion (Blaes 1987), have been studied. 

Interest in the PP instability waned in the 1990s when Balbus & Hawley (1991) pointed out that the Magneto-Rotational 
Instability (MRI), originally studied for magnetized Taylor-Couette flows (Velikhov 1959; Chandrasekhar 1961), can be im- 
portant for astrophysical accretion discs. Since the MRI is robust and requires only a weak magnetic field, the nonlinear 
development of MRI may lead to efficient angular momentum transport in accretion discs. Over the last two decades, numer- 
ous studies have been devoted to the MRI and related issues such as MHD turbulence in the disc (see, e.g., Balbus & Hawley 
1998 and Balbus 2003 for reviews; a sample of recent numerical studies include Hirose et al. 2009; Guan et al. 2009; Simon et 
al. 2009; Davis et al. 2010; Fromang 2010; Longaretti & Lesur 2010; Sorathia et al. 2010). 

Nevertheless, the question remains as to what happens to the original PP instability in an accretion torus when a finite 
magnetic field is present. After all, the tori produced in various astrophysical situations (e.g., binary mergers; Rezzolla et 
al. 2010; Montero et al. 2010) are expected to be highly magnetized. One might dismiss this question as purely academic since 
such a magnetic torus is likely MRI unstable and therefore turbulent. We note, however, that the usual MRI operates on 
perturbations with vertical structure (i.e., with finite vertical wavenumber k z ), while the PP instability operates on perturba- 
tions with k z — 0. That is, the PP instability pertains to the height-averaged behavior of the disc. Therefore one might expect 
that the PP instability will continue to operate even in the presence of MRI-induced turbulence. Furthermore, in connection 
with Galactic black-hole X-ray binaries, it has been suggested that accretion tori can support discrete, trapped oscillation 
modes, which might explain high-frequency quasi-periodic oscillations (e.g., Strohmayer 2001; Remillard & McClintock 2006) 
observed in a number of X-ray binary systems (e.g., Rezzolla et al. 2003; Lee et al. 2004; Schnittman & Rezzolla 2006; Blaes 
et al. 2007; Sramkova et al. 2007; Montero et al. 2007). Although it is currently not clear that pressure-supported tori provide 
a realistic model for the accretion flow around a black hole in any spectral state, structures resembling pressure-supported 
tori do appear to be present in some non-radiative global MHD simulations of MRI-driven turbulent accretion flows (e.g., 
Hawley & Balbus 2002; De Villiers et al. 2003; Machida et al. 2006). 

There have been a number of previous studies on global MHD instabilities in accretion flows. For example, Knobloch 
(1992), Kumar et al. (1994), Curry, Pudritz & Sutherland (1994) and Curry & Pudritz (1995) carried out global analysis for 
the axisymmetric modes with finite k z in differentially rotating flows threaded by vertical and/or azimuthal magnetic fields, 
thus establishing the robustness of MRI in these flows. Ogilvie & Pringle (1996) studied the non-axisymmetric instability of 
a cylindrical flow containing an azimuthal field, while Curry & Pudritz (1996) studied a similar flow containing a vertical 
field. Both studies focused on modes with finite vertical wavenumbers, which inevitably invite MRI. Although the effect of 
boundaries is emphasized, a somewhat arbitrary rigid boundary condition was adopted in these studies. As far we are aware, 
the behavior of the PP instability for finite tori with magnetic (both vertical and azimuthal) fields has not been clarified. 

In this paper, as part of our ongoing investigation of global oscillation modes and instabilities of rotating astrophysical 
flows (Tsang & Lai 2008, 2009a, b; Lai & Tsang 2009; Fu & Lai 2010), we carry out global stability analysis of magnetized 
accretion tori subjected to nonaxisymmetric perturbations. Since our main aim is to understand the effects of magnetic fields 
on the original PP instability, we focus on modes with no vertical structure (k z = 0) and we pay particular attention to the 
boundary conditions. As in many previous studies mentioned above, we model the torus by a cylindrical incompressible flow 
threaded by both vertical and toroidal magnetic fields. 

Our paper is organized as follows. In section 2, we describe the equilibrium models for our rotating magnetized flows. 
In section 3, the basic perturbation equations are presented. We derive the boundary conditions in section 4 and present our 
numerical calculations of the global instability in section 5. Final summary and discussion of our results are given in section 
6. 



2 EQUILIBRIUM MODELS 

As mentioned above, the PP instability operates in modes with no vertical structure (k z = 0). As such, the dynamics can 
be captured by height-averaged fluid equations. We consider a cylindrical shell (of finite width) of incompressible non-self- 
gravitating fluid, which is rotating differentially in the external gravitational field produced by a central compact object. We 
adopt the cylindrical coordinates (r, (f>, z) with the z-axis being the rotation axis. The cylindrical shell is assumed to be 
infinitely long in the ^-direction and threaded by magnetic fields. The fluid satisfies the ideal MHD equations: 

^ + (u-V)u = -ivn-V<&+-^(B- V)B, (1) 
dt p 4np 

^=Vx(nxB), (2) 

V • B = 0, (3) 
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Figure 1. The two magnetic field profiles adopted in the equilibrium torus model. 



V • u = 0. (4) 

Here, p is the constant fluid density, u the fluid velocity, B the magnetic field, and II = P + B 2 /8iv the total pressure with P 
being the gas pressure. The gravitational potential is $ = —GM/r, where M is the mass of the central object. The background 
flow velocity u — rQ(r)cj> and magnetic field B — B$(r)(j) + B z (r)z also depends only on r. For convenience, we assume the 
flow has a power-law rotation profile 

fi(r) oc r p . (5) 

The flow is confined between two boundaries (ri ^ r ^ T2), where the gas pressure vanishes (P\ T1 , r 2 = 0). Outside the 
fluid zone is a vacuum devoid of matter but maybe permeated with magnetic fields. In the equilibrium state, we assume that 
the magnetic field is continuous across the fluid boundaries so that there is no surface electric current at r = n, (however, 
we allow for surface current to develop when the fluid is perturbed). We will consider two models of magnetic field structure. 



2.1 Model (a) 

In this model, we assume that there is an external current running vertically at small radii much inside n, giving rise to 
B^(r) oc r" 1 in the inner region (r < ri). There is no azimuthal current in this region, so B z is constant. In the fluid zone, 
we adopt a power-law magnetic field profile 

B^(r)ocr', B,(r) oc r s , (6) 

which means that both the azimuthal and vertical current densities are also of power-law form. Outside the fluid zone (r > rz), 
there is no current. Hence, -B^(r) oc r _1 , B z (r) = const. The complete magnetic field profile is illustrated in the upper two 
panels of Fig. Q] 

Integrating the radial equilibrium equation 

1 -f=-^ + rtf-^, (7) 
p ar r z Aivpr 

gives the gas pressure profile 
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- - = 1 + - \*U ~ c > (8) 



P _ GM_ r z tt 2 _ 1 
7 ~~ + 2p + 2 ~~ 2' 

where 

vaz = Bz/^/Anp, v A <f> = B^/^jA-Kp (9) 
are the Alfven velocities and C is the integration constant. The location of gas pressure maximum is determined by 

~™ - " ^ " (i + ^ + ^ 2 - 0, (10) 
dr r z r r 

which defines a reference radius ro: 

^-=r n -s — - 1 + q — - . (11) 

r r 

Let C = XGM/ro with A being a constant and use Eq. (|11[) to substitute GM in Eq. (J5J, we can rewrite Eq. ([8| in the 
dimensionless form 

P - = \ - A + £^ + -Lo (-i + A - 1,-) + (1 + 3 )^ (-i + A - ir^) , (12) 
where 

VA4.0 = VA<t,{ro)/(r n ), v A zO = UA z (n))/(n)fio)- (13) 

Here and hereafter we will use units such that ro = fio = 1- Once we specify p, q, s, va</>o, vazO and A, we can determine the 
locations of the torus boundary by solving P = 0. However, there are several constraints on these parameters: 

(i) dP/dr = only guarantees the extremum of the P(r) profile. To ensure that we find a pressure maximum instead of 
minimum, we require d?P/dr 2 < at r = ro = 1, which implies 

2 P + 3 - s(l + 2s>Lo - (1 + g)(l + 2q)v A <t,o < 0. (14) 

This requirement reduces to p < —3/2 in the B — limit. 

(ii) Both sides of Eq. (|11[) need to be positive so that the gas pressure maximum exists. Thus 

1 - SVAzO - (1 + q)v A 4>Q > 0. (15) 

(iii) The maximum gas pressure P max must be positive. Thus, requiring the RHS of Eq. (|12[) to be positive at r — 1 gives 



A < 



1 v Az0 2v A(j>Q 



2p + 3 3 2 o 2 

7.VAzO — oVa^q 



2p + 2 2 



(16) 



provided that Eq. (|15|l is satisfied. 

Figure [2] illustrates some examples of Model (a). We specify the values of p, q, s, VA<f,o an d VazO, then by varying A, we 
obtain solutions for different torus thickness. For a given Ar/r2, both ri and change when VA<t>o changes as a result of 
magnetic support in the torus. In the hydro limit («a<^o = vazO = 0), n approaches 0.5 as Ar/r2 — > 1. This feature is shown 
analytically in Pringle & King (2007). For a finite field strength, we see that n — ► and f2 -> oo as Ar/r-2 —¥ 1. Note that 
for a relatively thin torus (Ar/r2 < 0.6), the differences of r2 and r\ between different field strengths are quite small. Since 
Btf, and B z have similar effects on the equilibrium structure (see Eq. [12]), these features also apply to models with a finite 
vertical field. The special case of s = and q = — 1 is worth mentioning: In this case, the magnetic field is force-free and has 
no effect on the equilibrium structure 0- 



2.2 Model (b) 



The magnetic field profile in this case is shown in the bottom two panels of Fig.[T] Compared with Model (a), the difference is 
that there is no vertical current at small radius. Therefore the azimuthal field in the inner region (r < n) is zero. In the fluid 
zone, we assume B z oc r and B^(r) oc r ~r\jr such that both azimuthal and vertical current densities are uniform. Following 
the same procedure as in section 2.1, we can derive the dimensionless expression for gas pressure profile: 



- - - - X 

p ~ r + 2p + 2 



r 2p+2 2 

+ V A zO 



\ 1 1 2 

A r 

r 2 



+ V A cj>0 



A-i 



1 + rf 



/ 2 

r(r 



2ri lnr) 



1-rf (1- 

Similarly, for a viable equilibrium model to exist, the model parameters must satisfy the following requirements: 



(17) 



1 This is why in Curry & Pudritz (1996) the one-to-one mapping between rijr\ and (r2 — ro)/r2 remains unchanged for different uniform 
vertical B field strengths. 
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Figure 2. Some examples of Model (a) with a pure toroidal magnetic field (B^ arm the fluid and B z = 0) and constant angular 
momentum distribution (Q, oc r~ 2 ). The x-axis is the relative thickness of the torus with Ar = T2 — r\ being the absolute width, and the 
y-axis shows the locations of the two boundaries. The different lines represent different values of va</>0 = v A<t>{ r o) / { r 0^o)i as indicated. 
The horizontal line indicates the location of gas pressure maximum (ro = 1). 



2p + 3 - 3v Az0 - 



3 + rf 

r 2\2 



(1 



«4<M) < 0, 



VAzO 



l+r 



12 ^ n 



A < 



- a 2 

v AzO 



T^r? A(j>0 



2p + 3 3 



2p + 2 



2 + rj a 

1 _ '2 V A4>0 



(18) 

(19) 
(20) 



3 MHD EQUATIONS FOR PERTURBATIONS 

Assuming that the Eulerian perturbation of any physical variable / is of the form Sf oc e lm 'l>~' tjt (with no dependance on z), 
the linearized perturbation equations are 

1 d , . . im . . 
- — (r<5u r ) + — <Sit = (21) 
r or r 

-iCj5u r - 2fiJ««, = -i^ + ^^B,. -Jit-SB* (22) 
p or 47rpr 27rpr 

-iw^ + ^tftir = ^m + -^(|- + iWtfB r + $^*fl* (23) 
v 20, pr 4np \dr r J v Anpr v 

. imB^ 1 dB z 

— iu)bu z = — -oB z + ; — 6B r (24) 

47rpr 47rp ar 

—iu}5B r = im — -Su r (25) 
r 

imBtt, d ( B<),\ dO foa\ 

—iluoBs — — oWrf, — r— — - ott r + r——oB r (26) 

r or \ r J ar 

-in imB^ dB z 

— ILOOB Z — 0U Z ; — ou r , (27) 

r ar 
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where Co = u — mQ, is the wave frequency in the co-rotating frame and the radial epicyclic frequency k is given by 

k 2 = — — (r 2 «) = 2(p + 2)fi 2 . (28) 
r dr 

Using An = Su + £ • Vu = d£/dt — —iuj^ + (u ■ V)£, we find that the Eulerian perturbation Su is related to the Lagrangian 
displacement vector £ by Su = — iui£ — rf2'£ r </> (prime donates radial derivative) and we can further combine Eqs. (|2HI - (|27ll 
into two equations for £ r (radial Lagrangian displacement) and 611/ p: 

^L =All ^ + Al2 — , (29) 
dr p 

d fsn\ „ ^ „ 5U 

dr V P J P 

where 

1 Co 2 — 2mCjQ, + m 2 wL 

An = - — ^, (31) 

r lo a — m z aj 

A U = ^, (32) 

A 2 i = w - m - 2rfi— + 2— — £ - 1 - 4— v (33) 

v dr \ dlnr ) v (cj 2 — m 2 uj A<t> ) 

2m Cba + mu 2 M 

A22 = t; , 2 , (34 J 

r u'- m ^a^> 

and cja^ = VA</>/r = B^l (ryf^Tp) is the toroidal Alfven frequency. Equations (|29[l and (|30|) are the same as Eqs. (119) and 
(120) (derived for a pure toroidal magnetic field) in Chandrasekhar (1961). 

Note that although we start with a mixed magnetic field B = B^ej) + B z z, the final Eqs. (|29l) and (|30l) do not contain 
B z . The reason is that the z-component only appears in Eqs. (|24|) and (1271) . which in fact can be decoupled from the other 
five perturbation equations. Indeed, using Eq. (|25[) to replace 8B r in Eq. (|24l) . and combining with Eq. (|27l) . we find 

(£ 2 - m 2 uj A4 ,) \SB Z + ^~j~$ u rj = 0- (35) 

In general, uj 2 — m 2 ^ 2 ^ 7^ 0. Comparing the above equation with Eq. (|27[) we have 8u z = 0. This is to be expected since the 
perturbed quantities are assumed to be independent of z. Also note that when the wave frequency ui is real, the coefficients 
-An, A21 and A22 are singular at 

Co 2 = m 2 Ld 2 A4> . (36) 
We shall call them the Magnetic Resonances (MRs). Obviously, they reduce to the corotation resonance when B^, = 0. 



4 BOUNDARY CONDITIONS 



In general, the boundary of any magnetized flow should satisfy the following conditions: 
[pu n ] = 0, 
[n ■ B] = 0, 



B 

P + pu 2 n + -L 
Sir 

B n Bt 



pUnllt 



4?r 



(37) 
(38) 

(39) 
(40) 



where n is a unit vector normal to the boundary surface, the subscript n and t denote the normal and tangential components, 
and the square bracket represents the difference in a quantity across the boundary (e.g., Schmidt 1979; Shu 1992). For the 
system we study in this paper, there is no radial background flow, we only need to consider Eqs. (I38|) - (|40[ ) with ii„ = 0. 
Obviously, with the magnetic field continuous across the boundaries and with no radial field component, our equilibrium 
models constructed in Sec. 2 already satisfy the boundary conditions. In the perturbed state, the boundary conditions read 



A [n ■ B] 
B n Bt 

47T 



= 0, 

= 0, 

= 0. 



(41) 
(42) 

(43) 
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Note that A[n-B] = [(An)-B] + [n-AB], and since both n and An are the same across the boundary[f], Eq. (|41|) simply becomes 
n ■ [AB] — 0, where we have used [B] = (as assumed in our model setup). Since (AB) r = (SB + £ ■ Vi?) r = 8B r — ^B^/r, 
and both £^ and B^ are continuous across the boundaries, we find 

[SB r ] = 0. (44) 
On the other hand, Eqs. (|41[) and (1421) combine to give 

[An] = 0. (45) 

The condition (|43|) is already satisfied because A [B„Bt] = [( AB n )Bt] + [B n ABt] = when the background field is continuous. 
We note that the perturbed magnetic field does not need to be continuous across the boundary. This means that there could 
be surface current induced by the perturbation. 

To implement the two boundary conditions (|44I) -(I45 | ). we need to calculate the perturbed magnetic field in the vacuum 
region (r < r\ and r > r%). This can be done by solving 

V x SB = 0, V • SB = 0. (46) 
Clearly, SB is a potential field SB = V* with * (also oc e ™<*-^*) satisfying 

V 2 * = 0. (47) 
The solution of Eq. ((47j) is 

* = Cir m + Car"™, (48) 
where Ci and C2 are integration constants. Requiring SB to be regular at r — > and r — > 00, we find 

SB r = C 1 mr m -\ 8B4, = Ci — r m , for r < n, (49) 



and 

SB r = -C 2 mr- m -\ SB^ = Ca^r"™, for r > r 2 . (50) 

The two constants C\ and Ci can be determined by using [SB r ] = 0. At r = ri, this implies Cimr" 1 " 1 = imB^r/r (see 
Eq. [25] for SB r inside the fluid zone). Thus, 

C 1 = iB^ r r- rn \ ri . (51) 
Similarly, 

C 2 = -iB^ r r m \ r2 . (52) 

Since the detailed realization of the boundary condition [All] = depends on the specific equilibrium model, we address the 
two models separately. 

4.1 Model (a) 

At the inner boundary r = ri, by using Eq. (|5ip . we find that the Lagrangian perturbation of the total pressure in the vacuum 
just inside n (i.e., r — r 1 -) is given by 

m MB ± _ m = m 

1 47T Anr Anr 

In the fluid just outside ri (i.e., r — ri+), we have 

" 2 ± _ 

' r 2 



An\ ri+ =sn + £ r 



2p+l _ ^40 _ J_ ^ _ su 2^ q g)^^) 



(54) 



Thus the condition [Ail] = at r = ri can be written as 



sn , 

P 



r 2p+l + 1 



H J (l - SVAzO - (1 + <?)VA0o) 



= 0, at r = r 1+ . (55) 



Note that in deriving the above equation, we have implicitly used [B^\ = and [£ r ] = 0. The same procedure yields the 
boundary condition at r = r 2 



2 Note that in the non-axisymmetric case, the perturbed surface normal vector An is not the same as r ; see Schmidt (1979) for a 
derivation of An. 
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Figure 3. The growth rate of Papaloizou-Pringle instability (in units of Qrj, the fluid rotation rate at the pressure maximum of the 
torus) as a function of the relative thickness of the torus for different values of m. The rotation profile is f2 oc r~ 2 . This figure is similar 
to Fig. 1 in Blaes & Glatzel (1986) and to Fig. 1 in Abramowicz et al. (1987). 



sn , 
p 



2p+l 



— (1 - sva z0 - (1 + q)v Atj , ) 



0, at r 



(56) 



4.2 Model (b) 

The derivation is similar to Model (a). In this case, the boundary conditions are 



5U , 

P 



sn 



2p+l 



2p+l 



~n 1 - V AzO 



1 



1+rj 2 

1 12 V A<j>0 



0, at r = r ± + , 



(1 



1 I, 2 

1 - v AzO - 



1 + rj 2 
1 Z2 V M0 



0, at r 



(57) 
(58) 



5 NUMERICAL RESULTS 

We employ the standard shooting method (Press et al. 1992) to solve Eqs. (|29l) and (|30|l subjected to the boundary conditions 
(I55|l - (|56[ ) [Model (a)] or (I57|) - (|58|) [Model (b)] to obtain the eigenvalue uj = uj r + iuji. For most of our analysis, we set the 
rotation index p — —2 (i.e., Q. oc r -2 ) such that our results can be directly compared with the original Papaloizou-Pringle 
instability. 

Before discussing our results for finite magnetic fields, we briefly review the main features of the classical (hydrodynamical) 
PP instability. As seen in Fig. [3] the instability growth rate increases with increasing torus thickness for small Ar/r2 but 
terminates at some finite thickness. As m increases, the termination point shifts to smaller Ar/r2, although the peak growth 
rate remains approximately the same. This means that the PP instability only exists for relatively thin tori, as shown by Blaes 
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Ar/r 2 

Figure 4. The instability growth rate as a function of the relative thickness of the torus for a range of toroidal magnetic field strengths. 
Different lines represent different Vaso with the solid line denoting the hydrodynamic case. The upper and bottom panels are for the 
m = 1 and m = 2 modes, respectively. 

& Glatzel (1986) and by Abramowicz et al. (1987). The former also provides an approximate analytical expression for the 
limiting maximum growth rate as m — > oo. These features can be understood from the fact that the PP instability arises from 
the interaction of the surface gravity waves at the torus boundary. For a thin torus, the velocity shear across the corotation 
point is small, and there would not be enough shear rotational energy available to drive the growth. When the torus thickness 
is too large, the wave amplitudes at the corotation radius (where the waves exchange angular momentum) is too small to allow 
for adequate interactions. Thus, only for the "intermediate" torus thickness, with mAr/r ^ 1, will the instability operates. 
This explains why for a larger m, the PP instability terminates at a smaller torus thickness. 

5.1 Model (a): Pure Toroidal Field Configuration 

In this section, we present the numerical results for an accretion torus with a pure power-law profile toroidal magnetic field. 
We choose the power-law index q = 1 so that the vertical current density is uniform as in Model (b). 

Figure 2] shows the growth rate uJi as a function of thickness Ar/r2 for different va^o- The two panels share similar 
characteristics: (i) For relatively weak B field (va$q % 0.1), the instability resembles the 5 = limit in that it always starts 
from infinitely small thickness and terminates beyond a certain Ar/r?; (ii) For stronger B fields, the instability starts beyond 
certain finite Ar jr-i and then extends all the way to very large thickness (although as Ar/r2 approaches unity, the growth 
rate becomes increasingly small). As va</iO increases, the critical thickness for the onset of instability also increases. 

Figure [5] maps out the unstable zone in the thickness - magnetic field strength parameter space. It shows the similar 
feature as as Fig. [4] We can see that the unstable region is mainly located at the lower-left (thin torus with weak B field) and 
the upper-right (thicker torus with strong B field) corners of the parameter space. 

In Fig. [6] we present our numerical results in a different way. We fix the dimensionless thickness Ar fr-i and plot the 
growth rate as a function of magnetic field strength. For a thin torus, we find that as va<po increases, the growth rate first 
goes up slightly compared to the B = case, then decreases and becomes completely suppressed when the magnetic field is 
sufficiently strong (vA<t> comparable to rotation velocity). For a thick torus, the instability can survive for a relatively stronger 
B field and then vanishes beyond a certain «a</>o- 

To probe the underlying physics of how magnetic fields affect the PP instability, we show in Fig.[7]the locations of several 
special points in the fluid: The corotation radius r c is where the wave pattern corotates with with the background flow, i.e., 
tb r = ui r — rafi = 0. The inner/outer magnetic resonances (IMR/OMR) are defined by [see Eq. (|36[) ] 
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Figure 5. The instability region in the parameter space defined by the dimensionless torus thickness Ar/r2 and the toroidal magnetic 
field strength. The other parameters are fixed to m = 2, f! oc r~ 2 and B,p oc r. The dotted area denotes the region where a growing mode 
can be found. 




Figure 6. The instability growth rate as a function of = v Acf>( r o) / { r 0^o) f° r tori with different thickness Ar/r2. The upper and 

bottom panels depict the cases with m = 1 and m = 2, respectively. 
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Figure 7. Some special radii for the m = 1 overstable mode in a torus as a function of the dimensionless magnetic field strength VAcf>o- 
The upper and bottom panels correspond to a thin Ar/r2 = 0.2) and thick (Ar/ri = 0.7) torus, respectively. The two solid lines show 
the inner and outer torus boundaries (ri and r%)- The dotted line represents the corotation radius, while the shot-dashed and long-dashed 
lines show the inner and outer magnetic resonances (IMR and OMR), respectively [see Eq. {59}]. The thin solid line is the mode growth 
rate with the scale shown on the right. The torus rotation and magnetic field profiles are f! oc r~ 2 and oc r. 



UJr = OJr — mQ, — ±mUJAcj>- (59) 

At the IMR, the wave is trailing the background flow but corotates with the azimuthal Alfven wave traveling in the counter- 
rotational direction (viewed in the corotating frame), while at the OMR, the wave is leading the background flow and corotates 
with the Alfven wave in the rotational direction. Recall that for PP instability to operate in the B = limit, it is essential 
that the corotation radius lies in between torus boundaries (i.e., r\ < r c < ri). Now, with the inclusion of the magnetic 
field, we see from Fig. [7] that as VA(po increases, both r c and roMR shift beyond the outer boundary of the torus. The IMR 
radius, riMR, however, always stays inside the fluid. This suggests that in a magnetic torus, the IMR plays a similar role as 
the corotation resonance does in a non-magnetic torus. 



5.2 Other Magnetic Field Configurations 

5.2.1 Model (a): Mixed magnetic field 

Although the vertical magnetic field B z does not enter into the perturbation equations, the presence of a finite B z can affect 
the mode growth rate through the boundary conditions. Figure [8] shows some results for the accretion tori with a mixture of 
vertical and toroidal magnetic fields. We take the power-law index of B z to be s = 1 so that the azimuthal electric current 
density is constant. In Fig. [8] we fix the toroidal field and plot the mode growth rate as a function of the vertical field strength. 
We find that the effect of finite vazO is small (note the scale of the y-axis). For vaz = 0, the results agree with what is shown 
in Fig. [6]d. 



5.2.2 Model (b): Pure toroidal field 

In this case, since the background toroidal magnetic field has a profile that depends on the inner boundary radius ri, to 
solve for the equilibrium structure and the global mode, we must first specify ri. Once we fix n, we can easily solve for the 
other boundary radius ri. In Fig. [9l we show in the upper panel the mode growth rate as a function of va4>q- The result is 
qualitatively similar to Fig. [6p. The bottom panel shows that for a fixed inner boundary radius n, the thickness does not 
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Figure 8. The instability growth rate as a function of vazO f° r a torus with relative thickness Ar/r2 = 0.2. Different lines represent 
different values of t>ytc/>0- The other parameters are fixed to m = 2, f! oc r~ 2 , oc r and B z oc r in the fluid. 
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Figure 9. The instability growth rate as a function of vam for accretion tori described by Model (b) with a pure toroidal magnetic 
field. The solid and dotted lines correspond to two different inner disc boundary radii. The bottom pane shows the corresponding torus 
relative thickness. 
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Figure 10. The instability growth rate as a function of vazO f° r a thin torus (with r\ = 0.9 and Ar/r2 — 0.2). The different lines 
represent different values of VAtf,a- Unlike Fig. [8] here the magnetic field profile is described by Model (b). 



change appreciably as VA<f,o varies. So the two values of r\ we choose adequately depict the thin and thick tori, respectively. 
Again, we see that for a thin torus, the original Papaloizou-Pringle instability is suppressed by the toroidal field, while for a 
thick torus, the instability can survive for larger field strengths. 

5.2.3 Model (b): Mixed magnetic field 

In Fig. I lOt we show the mode growth rate as a function of VazO for different fixed values of VA<j>o- Similar to the case shown 
in Fig. [5] we see that the vertical magnetic field has a small effect on the stability property of a magnetized accretion torus. 
Again, this is understandable given that B z does not enter into the differential equations for the perturbations, but only 
affects the modes through boundary conditions. 



6 DISCUSSION 

In this paper, we have studied the effect of magnetic fields on the global non-axisymmetric instability (the PP instability) in 
accretion tori. For simplicity, we assume that both the perturbation and the background flow variables have no z-dependance 
(thus our tori are essentially 2D cylinders). We have explored various possible magnetic configurations in the torus. Although 
the detailed property of the instability is model-dependent, Figs. |4][6] illustrate our general findings: (i) For thin tori (with the 
dimensionless thickness Ar/r2 ^ 0.2, where V2 is the outer torus radius), the instability exists for zero and weak magnetic 
fields, but is suppressed when the toroidal field becomes sufficiently strong (with the corresponding Alfven speed va</> <; 0.2rfi 
measured at gas pressure maximum); (ii) For thicker tori (Ar/r2 <; 0.4), the PP instability does not operate for zero and weak 
magnetic fields, but becomes active when the field is sufficiently strong (va$ ii 0.2rQ measured at gas pressure maximum). A 
vertical magnetic field also influences the PP instability, but its effect is generally smaller that that of the toroidal field. 

It is difficult to precisely pin down the physical origin of the magnetic field effect on the PP instability. For example, with 
a finite toroidal magnetic field, we find that the corotation resonance radius may lie outside the torus body, and yet the torus 
is still unstable. On the other hand, the inner magnetic resonance radius, where ui r — mfi = —mu)A<t> [see Eq. (|59[)]. always 
lie inside the fluid body. Thus we suspect that in a magnetic torus, the inner magnetic resonance plays a similar role as the 
corotation resonance does in a non-magnetic torus. 

We note that the PP instability (or its magnetic generalization) involves wave modes that do not have vertical structure 
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(i.e. k z = 0). Thus it is distinctly different from the usual MRI (Balbus & Hawley 1998). Our finding about the instability of 
think tori with strong magnetic fields is particularly interesting: Since the MRI can be suppressed when the magnetic field is 
too strong, our results suggest that magnetized tori may be subject to the instability even when it is stable against the usual 
MRI. 
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